arXiv: 1501.07804v6 [nlin.CD] 29Jun2016 


Lagrangian transport through surfaces in 
volume-preserving flows 

Daniel Karrasch*^ 

ETH Zurich, Institute for Mechanical Systems 
Leonhardstrasse 21, 8092 Zurich, Switzerland. 

June 30, 2016 


Advective transport of scalar quantities through surfaces is of fundamental 
importance in many scientific applications. From the Eulerian perspective of 
the surface it can be quantified by the well-known integral of the flux den¬ 
sity. The recent development of highly accurate semi-Lagrangian methods for 
solving scalar conservation laws and of Lagrangian approaches to coherent 
structures in turbulent (geophysical) fluid flows necessitate a new approach 
to transport from the (Lagrangian) material perspective. We present a La¬ 
grangian framework for calculating transport of conserved quantities through 
a given surface in n-dimensional, fully aperiodic, volume-preserving flows. 
Our approach does not involve any dynamical assumptions on the surface or 
its boundary. 


1. Introduction 

The transfer of a quantity along the motion of some carrying fluid, or Lagrangian trans¬ 
port for short, is of fundamental importance to a broad variety of scientific fields and 
applications. The latter include (geophysical) fluid dynamics [30, 37, 34, 20], chemical 
kinetics [ ], fluid engineering [ , ] and plasma confinement [ ]. Existing methods 

for computing transport (i) have mostly been developed under certain assumptions on 
temporal behavior (steady or periodic time-dependence) [23, 24, 22, 3 |; spatial location 
(regions related to invariant manifolds such as lobe dynamics and dividing surfaces in 
transition state theory) [23, 24, 22, 33, 25, 1]; state space dimension (2D) [23, 33, 25, ]; or 
(ii) restrict to a perturbation setting | ]. Recently, the problem of quantifying finite-time 
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transport in aperiodic flows between distinct, arbitrary flow regions has been consid¬ 
ered by Mosovsky et al. \ ]. They present a framework for describing and computing 

finite-time transport in n-dimensional (chaotic) volume-preserving flows, which relies on 
the reduced dynamics of an (n — 2)-dimensional ‘minimal set’ of fundamental trajecto¬ 
ries. In this paper, we present a Lagrangian approach to the complementary problem 
of computing transport through a codimension-one surface over a finite-time interval in 
volume-preserving flows. This cannot be reformulated in the setting of [ s] ; since (i) 
initial and final positions of surface-crossing particles are a priori unknown, and (ii) par¬ 
ticles may cross the surface several times, possibly in opposite directions, leading to a 
net number of surface crossings different from one, in particular including zero. 

The problem of computing the flux through a surface in general flows admits a well- 
known solution in terms of an integral of the flux density over the surface and the time 
interval of interest, cf. the left-hand side of Eq. (1). We view this approach as Eulerian, in 
that it involves instantaneous information (flux density) at fixed locations in spacetime. 
Recently, two lines of research emerged that inevitably require a Lagrangian approach to 
the flux calculation. 

The first is concerned with the numerical solution of advection equations (or, in the 
absence of sources, conservation laws) for conserved quantities by semi-Lagrangian meth¬ 
ods, which enjoy geometric flexibility and the absence of Eulerian stability constraints 
[38, 31 1. Roughly speaking, the term semi-Lagrangian refers to methods which evolve 
material densities on spatially fixed test volumes based on short-term material advection 
steps. 

The second is concerned with transport by coherent vortices in oceanic flows, more 
precisely with the determination of the relevance of coherent transport [6, 42, 17], an 
open problem in physical oceanography and climate science. Coherent structures have 
long been studied in fluid dynamics [7], typically from an Eulerian point of view, i.e., by 
defining coherent structures as subsets of spacetime , usually as (sub-)level sets of scalar 
fields on spacetime [26, 32, 6, ]. This approach, however, yields coherent structures 

with a priori unclear relation to actual fluid motion and hence ambiguous role in co¬ 
herent transport; see transfer-operator related approaches which seek Eulerian coherent 
sets with minimal flux through the boundary under advection with small-scale diffusion 
superimposed, cf., for instance, [9]. For these reasons, Lagrangian approaches to coherent 
structures have been developed over the last few years. These seek coherent structures 
as subsets of the material evolving under the flow [14, 15, 18, 13, 10, 11, 21, 29], building 
on a wide variety of mathematical principles. The location of a material structure in 
spacetime is fully determined by its motion. 

This implicit spatial definition of material structures poses a severe challenge to deter¬ 
mining their contribution to transport through a surface in an Eulerian manner. First, 
for each point on the extended surface in spacetime knowledge is required whether it is 
occupied by a particle originating from the material structure of interest. Second, the 
subset of points on the surface which are occupied by material particles of interest may be 
very complicated, cf. Fig. 1, especially in turbulent flows. As a consequence, an Eulerian 
integration of the flux density restricted to the intersection of the extended surface with 
the path of the material set of interest is practically infeasible. 
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Figure 1: Views from opposite directions on the evolution of an initial ellipse over time, 
and its intersection with a surface in extended state space. Computing the flux 
contribution by the material points originating from that ellipse in an Eulerian 
framework corresponds to computing the flux integral over the subset of the 
section which is contained in the interior of the cylindrical structure. 


Conversely, from the material point of view, it is intuitive that the number of sur¬ 
face crossings of each individual Lagrangian particle is relevant to determine its flux 
contribution. First steps towards a classification of Lagrangian particles with respect 
to net number of curve crossings have been proposed for two-dimensional flows under 
additional technical assumptions by Zhang [10, ]. In this paper, we solve the following 

more general donating region problem [39, Def. 1.2], 

Problem 1 For a given regular, divergence-free and time-dependent vector field u{t,x) = 
uflx), consider a conserved quantity f(t,x) = ft(x), which satisfies the scalar conserva¬ 
tion law in Eulerian form 1 

d t f + u t ■v / = o. 

Let C be a compact, connected, embedded codimension-one surface in (configuration space) 
M n ; and T = [0, r] be a compact time interval. The problem is to find pairwise disjoint 
sets T>f; C M n , indexed by k G Z, of Lagrangian particles at time t = 0, such that 

[ ft(x) ■ u t (x) • n(x) dxdi = Y'b [ ft\ t=0 (p ) dp, (1) 

JCxT k£Z Jv k 

where n is the unit normal vector field to C characterizing the direction of positive flux. 

In fact, we may generalize Problem 1 and drop the assumption that C be stationary over 
T. Thus, we solve the analogous problem for smoothly moving surfaces TL = {J te q-Ct, 
i.e., TL is everywhere transversal to time fibers {t} x M n . 

Problem 2 Let u(t,x ) = uflx) and f(t,x ) = ft(x) be as in Problem 1 , and uflx) = 
(ut{x) ) extended velocity field. Let Tl be a codimension-one surface in extended 

lr The Lagrangian form is = 0 with flf the material derivative. Physically, this means that the scalar 
/ does not change along particle motions. 
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configuration space T X M n ; T = [0,r], as defined above. The problem is to find pairwise 
disjoint sets T >C M n , indexed by k € Z, of Lagrangian particles at time t = 0, such that 

[ fi(x) • u t (x) • ht(x)d(t,x) = V k • f ft\ t=0 (p)dp, (2) 

where nfix) is the unit normal vector field to TL at (t,x). 

For ease of presentation, we will first focus on the stationary case as in Problem 1, 
before we discuss the general case as in Problem 2 in Section 5. 

2. Illustrative discussion of the main result 

The aim of this section is to discuss Eq. (1) in anticipation of mathematical construc¬ 
tions presented in Sections Section 3 and Section 4. In fact, those constructions emerge 
naturally when taking a differential topological view on the change from Eulerian to La¬ 
grangian coordinates. We recall and discuss concepts related to that coordinate change 
in Section 3. 

First, the left-hand side of Eq. (1) corresponds to the classical flux integral over the 
surface C over time T. The integrand, also called the flux density, corresponds to the 
normal component of the velocity field weighted with the material density /. Here, the 
flux is viewed from the perspective of the section, which is why we refer to the left-hand 
side of Eq. (1) as the Eulerian flux integral. 

In contrast, the right-hand side of Eq. (1) is an integral, or a sum of integrals, solely 
over material points or particles. As our later analysis reveals we may decompose the 
whole material domain (up to a measure zero set) into disjoint material subsets, each of 
which contains particles with a certain number of net crossings across the surface C over 
the time interval T. These material subsets are denoted by D*,, where k refers to the net 
number of transversal crossings, and have been coined donating regions of fluxing index 
k by Zhang [4 ]. The contribution to transport by each donating region corresponds to 
its respective mass (or volume in a homogeneous fluid) multiplied by its corresponding 
fluxing index. The overall transport through C is then given by the sum of transport 
contributions of all donating regions. In summary, the right-hand side of Eq. (1) views 
the flux from the perspective of the mass-carrying particles, which is why we refer to it 
as the Lagrangian flux integral. 

For illustration, consider an array of vortices, given by the stream function 

H(x,y) = —A sin (nx) sin (ny ), 

uniform spiraling forcing given by 
= (t sin ( 7 rt ), t cos (7rt)) T . 

_ f d y H(x,y) — t sin (irt)\ 

\-d x H(x, y) - t cos ( 7 rt )) ' 


subject to an aperiodic, spatially 

F(t) 

The velocity held is hence 

u (t,x,y ) 
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Figure 2: Donating regions divided in cells with fluxing index +1 (vertically hatched), — 1 
(diagonally hatched) and 0 (white). The dashed curve (red) is a sample trajec¬ 
tory (path line). The ellipse (dark green) corresponds to the set of Lagrangian 
particles whose flow evolution is shown in Fig. 1. 


We set the section to C = {x = 0.75, —0.2 >y> 1.2}, the time interval to T = [0, 2.5], 
and, for simplicity, the material density / = 1. Fig. 2 shows the regions of particles with 
fluxing index +1, 0, and —1; here, white corresponds to fluxing index 0. For details on 
how to construct the donating regions, see Section 4.4. There are no particles with other 
fluxing indices, and some particles attain a fluxing index +1 by crossing twice in positive 
and once in negative direction, such as the sample trajectory shown in red. 

For a numerical validation of Eq. (1), we compute its left-hand side by numerical 
quadrature, which yields a transport value of —0.142. For the right-hand side, we subtract 
the area of D_i from that of T> i, and obtain a transport value of —0.141, with relative 
error below 1%. For further improvements of area calculations via spline-approximation 
of polygon boundaries see [38]. Note that, in practice, the numerical integration of the 
Eulerian integral becomes more challenging in case the section C is curved and requires a 
parametrization. In contrast, we will see that a curved section does not add any difficulty 
in the Lagrangian framework. This geometric flexibility has been one major motivation 
to develop semi-Lagrangian numerical schemes for solving conservation laws, cf. [38, ]. 

For our main intended application, the calculation of the contribution to transport 
by a certain material subset A —such as the ellipse in Figs. 1 and 2 (dark green, cen¬ 
tered at (0.3,0.5) with semi-major axis 0.6 in y-direction, and semi-minor axis 0.4 in 
x-direction)—reduces to subtracting the area of A n P_i from 4 n Di, which yields 
-0.1061. 

At first sight, this may seem paradoxical, since the entire ellipse is launched from one 
side of the section. Fig. 1 now indicates that at the end of the time interval, parts of 
the material ended up on the other side of the section. Going from the first to the latter 
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side corresponds to positive flux, and one would expect an overall positive transport 
contribution. However, a closer inspection of Fig. 1 reveals that most of the material 
which is advected to the other side has actually not crossed the section, but has flown 
around it. Moreover, much of the material which has flown around has actually crossed 
the section in negative direction, thus yielding an overall negative flux contribution of 
the ellipse. 

The rest of this paper is devoted to the derivation of Eqs. (1) and (2). Most of the 
involved abstract concepts turn out to have a nice physical interpretation. 

3. Eulerian and Lagrangian coordinates 

While Eulerian coordinates x are assigned to spatial points in a fixed frame of reference, 
Lagrangian coordinates p label material points and are usually taken as the Eulerian 
coordinates at some initial time, say, t = 0. The motion of material points is described 
by the flow ip, a mapping between initial positions p at time t = 0 and current positions 
x at time t. i.e., Pq(p) = x. Thus, the flow map cpg can be interpreted as a change from 
Lagrangian to Eulerian coordinates. The inverse flow map iff = (pq) corresponds then 
to the change from Eulerian to Lagrangian coordinates. 

In fluid dynamics, there are two important characteristic curves associated with the 
flow [ ], which we re-interpret in terms of Eulerian and Lagrangian coordinates. 

A path line through p is the time-curve of a fixed Lagrangian particle p in Eulerian 
coordinates, i.e., t <—> tpy(jp). In other words, the path line is a collection of Eulerian 
positions that the Lagrangian particle p will occupy at some time. Its time derivative, 
expressed in Eulerian coordinates, gives rise to the velocity field, 

u t {x) ■■= u(t,x ) = d t p o(^(x)). (3) 

By construction, the path line through p is the solution of the initial value problem 
x = Ut(x), x(0) = p. 

A streak line through x is the time-curve of a fixed Eulerian location x in Lagrangian 
coordinates, i.e., t i —> p^fix). In other words, the streak line is a collection of material 
points that will occupy the Eulerian position x at some time, see Fig. 3. Our definition 
of streak lines suits well the currently intended purpose and is consistent with Zhang’s 
use [ ]. It is more common and actually more intuitive, however, to view a streak line as 

the collection of material points that have passed the Eulerian position x at some time. 
In this context, a streak line can be imagined as an instantaneous curve of Lagrangian 
markers, injected in the past at x and passively advected by the flow, see [ ]. The time 
derivative of the streak line, expressed in Lagrangian coordinates, gives rise to the streak 
vector field 

w t (p) := w(t,p) = dtPt(<Po(p))- 

By construction, the streak line through x is the solution of the initial value problem 
p = w t (p), p(0) = x. 
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Figure 3: Streak line and streak vector field: the streak line through xq (red) is con¬ 
structed by subsequently advecting xq from time t backwards to time 0. Con¬ 
versely, when a particle from the streak line, with curve parameter s (colorbar), 
is launched at time 0 and advected s time units, the particle reaches the spatial 
position .To (black). The streak vector field corresponds to the velocity of the 
streak lines. The bottom x-y-plane is to be considered as Lagrangian coor¬ 
dinates. Correspondingly, the streak vector field is a material vector field, in 
contrast to the classic velocity field, which is a spatial vector field. 


To derive an algebraic relation between the velocity field u and the streak vector 
field w , we compute the velocity dtPt(x) along the streak line t K > Pt{x) through x. 
Differentiation of the constant function 1 H > p^ (p® (*)) = x with the shorthand notation 
V = <Pt ( x ) yields 

0 = d t Po (p) + dpQ (p) d t p° t (s). 

Using Eq. (3) and the invertibility of the linearized flow map with dpQ (p) -1 = dp^ (x). 
we obtain 


w t {p) = - (d^ (p)) 1 u (t, Po (p)) . (4) 

Equivalently, in global terms we have 

w t = - (<£?)* u t, u t = -(p t 0 ) if w ti 

where the subindex * denotes pushforward of vector fields by p® and p q, respectively. 
Eq. (4) is an alternative to the formula originally derived by Weinkauf &; Theisel [36]. 

Finally, we recall that the flow p is volume-preserving on W 1 if and only if the velocity 
field u (or, equivalently, w) is divergence-free [2]. 

4. Transport through surfaces 

In this section, we solve Problem 1 by analyzing the change from Eulerian to Lagrangian 
coordinates from the differential topology viewpoint [27, 16]. 
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Figure 4: Eulerian coordinates: Flux density on the extended section 7~L, bounded by 
blue, yellow, purple and brown lines. Positive flux is directed in positive re¬ 
direction, and the flux direction through % is indicated by +/—. Zero-level 
curves of flux are shown in black. Also shown is the extended pathline from 
Fig. 2 (red, dashed when behind the section) with three transversal crossings, 
two positive and one negative. 

4.1. Setting 

We decompose the change from Eulerian to Lagrangian coordinates (evaluated along the 
extended section T xC) into two steps. First, we map Eulerian spacetime points to their 
respective initial position at time t = 0, while keeping them on the same time slice, i.e., 

<F: T x R n -> T x M n , (t,x) i-> (t,p?(s)) = (t,p). 

The transformation $ maps the (extended) section 7L :=T xC, see Fig. 4, diffeomorphi- 
cally to the streak surface S := <!>(%), see Fig. 5. 

In extended Eulerian coordinates, trajectories (or, extended path lines) take the form 
(t,</?Q(p)); in extended Lagrangian coordinates, trajectories are simply vertical lines, i.e., 
t i— > ( t,p ), since they track the same particle p, see Fig. 5. Intersections of path lines 
with 7~L correspond one-to-one with intersections of S with vertical lines. 

In a second step, we project spacetime points to their spatial coordinates by the canon¬ 
ical projection II : T xM n —> M n and introduce V := 11(5). We view the image of II as the 
initial time slice of the extended state space and therefore parametrized by Lagrangian 
coordinates. By construction, we have for T := II o <F the identity T(t,x) = p^x), and 
it is exactly the particles from T> that cross C within T one way or another, possibly 
multiple times. 

For later reference, we emphasize the following observation. Let be the unit time-like 
vector held. Then we have 


d(/>od'F(t, x)e t (x) = -u t (x). 


(5a) 
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Figure 5: Lagrangian coordinates: Streak surface S and trajectory (red), which crosses 
S twice in positive (blue balls) and once in negative (red ball) direction. The 
boundary curves correspond to the boundary curves in Fig. 4 with the re¬ 
spective colors. Their projections onto the t = 0 time slice are shown in the 
corresponding colors, cf., Fig. 2. 

On the other hand, for any space-like vector held Vt such as, for instance, the velocity 
held u, we have 


d^gd ^(t,x)v t (x) = v t (x), (5b) 

since in this case dT acts like 6.$. 

4.2. Counting net crossings — the degree 

We study the differentiable map \L as a map between the two n-dimensional manifolds 
WcTx M n and T> C M n . To this end, we recall notions from differential topology and 
interpret them in our setting. 

First, regular points (f, x) £ 7i are those for which the differential d\l/(t, x) is invertible. 
This is exactly the case, when the tangent space at <F(f, x) = to S does not 

contain the vertical time direction. Non-regular points (t, x) of T are referred to as 
critical points, which are characterized by non-transversal crossings, i.e., instantaneous 
stagnation points ut(x) = 0, or tangential crossing/rebound. 

According to Sard’s theorem, the set of critical values, i.e., images of critical points, has 
measure zero in V, even though the set of critical points may be large in 7i. Consistently, 
critical points do not contribute to the Eulerian flux integral. 

In the Eulerian setting, positive/negative crossings of a Lagrangian particle p through 
C correspond to ±«t(a:)-n(x) > 0 at corresponding intersections (t, x ). In the Lagrangian 
setting, positive/negative crossings are equivalent to crossings of the vertical line over 
p through S from below/above with respect to the orientation on S, see Fig. 5. This 
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is formalized by TdetdT > 0, which geometrically means that d\I f(t,x) is orientation¬ 
preserving or orientation-reversing at the respective crossings of p. This is well-defined 
at all crossings exactly for particles p which are regular values of T. 

The difference of positive and negative crossings, or, in other words, the net number 
of crossings, is given by the degree of T at p £D, 

deg(T,p) := E sign (det (dT(t, x ))). 

( t,x)£'&- 1 (p)r\'H 

The degree function is locally constant, and can be used to define a partition of the set of 
particles, up to the measure-zero set of critical values. Specifically, we define T>k as the 
set of regular values p of T such that deg(T,p) = k for k £ Z. Following [41], we refer 
to the T>k s as donating regions of fluxing index k. By construction, each T>k contains 
the particles with only transversal crossings of which there are k net. In particular, all 
material points that do not cross C within T are contained in Dq. 


4.3. Area formula 

In the homogeneous case / = 1, we have the following identity for the degree of T 
[12, Section 3.1.5, Thm. 6], which is a consequence of the area formula ([8] and [19, 
Thm. 5.3.7]): 



detd\k(t, x) dx d t 


deg(T,p) dp = ^2 &vol(Pfc). 


The area formula is a generalization of the change of variables (in integral) formula to non- 
injective maps such as T. Clearly, the degree measures the non-injectivity, and is equal 
to ±1 globally for orientation-preserving and -reversing diffeomorphisms, respectively. 

Inhomogeneous conserved quantities are constant along trajectories in volume-preserving 
flows. Therefore, we have ft(x) = /o(4/(t, x)), and we conclude with the general area 
formula: 



ft{x) det dT(f, x) dx dt 



,/o( v 3/(t, x)) det dT(f, x) dx d t 



/o(p)dp. 


To prove Eq. (1) and thereby solve Problem 1, it remains to show 


detdT(t, x) = u t (x) • n(x), 


( 6 ) 


which is done in Appendix Section A, taking advantage of the relation (4) between the 
streak vector field and the velocity field. 


4.4. The two-dimensional case 

Our results for volume-preserving flows on M n simplify in the two-dimensional case con¬ 
sidered by Zhang [40, 41] as follows. It is well-known that, in our context, the degree 
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of 'I' at some regular value p G T> C M 2 equals the winding number of p with respect to 
the closed curve around p [5, Section 6.6]. The winding number counts the net 

number of turns of 'fy(d'H) around p under one counter-clockwise passage through dTi , 
see Fig. 2. 

In the 2D case, we find the donating regions 'Df ; of fluxing index k as follows. The 
boundary dT~L of the extended section % is a closed curve, just as its image under T. In 
Figs. 5 and 2 its image corresponds to the concatenation of the section (blue), one streak¬ 
line (say, yellow), the backward image of the section (purple) and the other streakline 
(brown). This closed curve gives rise to possibly several connected components, simple 
loops, through its self-intersections. Now, the winding number is constant on the interior 
of each simple loop, and it suffices to compute the winding number for some contained 
sample point, which is a standard task in computational geometry. As anticipated in 
Section 2, in the case of constant density the Lagrangian flux calculation reduces to 
the computation of the enclosed area of each component, multiplication by its winding 
number and final summation over the loops. 

In summary, we have generalized the main result of [ ], in which only the issue of 

counting net crossings was treated under additional assumptions from a different method¬ 
ological perspective. Note also that for the isolated counting aspect, our characterization 
in terms of deg (\L), or winding number in 2D, is also valid in the case of a compressible 
velocity held u. 

5. Moving surfaces 

To find a Lagrangian flux integral in the moving surface case, Problem 2, we may proceed 
as in Section 4 with % = T x C replaced by % = [J f Ct as the domain of the maps <F and 
\L. All constructions and arguments work just as well as before. It remains to find an 
expression analogous to Eq. (6) for detdT(f, x): 

detdT(f,x) = iit{x) • h t (x), 
which is shown in Appendix Section B. 

6. Conclusion 

In this work we have devised a Lagrangian approach to transport through (codimension- 
one) surfaces in general n-dirnensional, unsteady, volume-preserving hows. Studying 
the change from Lagrangian to Eulerian coordinates (and vice versa) alone, we have (i) 
discovered a striking analogy between path lines and streak lines, (ii) derived an algebraic 
relation between their associated velocity helds, (iii) found a natural way of determining 
the net number of surface crossings for individual Lagrangian particles up to a measure 
zero set, and, finally, (iv) transformed the Eulerian hux integral into a Lagrangian one. 
Thus, we have solved the donating region problem [39] for volume-preserving hows in 
arbitrary Unite dimension and smoothly moving surfaces. An obvious extension of the 
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methodology is to handle non-volume-preserving flows, and work is in progress in this 
direction. 

Besides the theoretical insights, the Lagrangian approach to transport through surfaces 
is of major importance in at least two relevant fluid mechanical applications. First, it 
provides the missing theoretical foundation of a family of highly accurate semi-Lagrangian 
finite volume and interface tracking methods [38, 39]. Second, it facilitates the efficient 
computation of transport by Lagrangian coherent vortices in large-scale oceanic flows 
[18, 17]. 

A. Proof of det = u • n 

We show that det d\L(t, x) = u t (x) • n(x). To this end, fix a regular point (t,x) E H. 
Then T acts diffeomorphically between an open neighborhood U of (t, x) in T~L and its 
image V := in V. We may choose local coordinates on U such that in (t, x ) we have 
an orthonormal basis, i.e., e* E TfT. span {e 2 ,..., e n } = T X C and e.\ E Tjr'H. with 

n 

f\ e.i = e 1 A ... A e n , e t A e 2 A ... A e n , and e t A e\ A ... A e n , 
i=l 

unit volume parallelepipeds. To avoid a discussion on delicate orientation issues, the fol¬ 
lowing calculations involving determinants are to be read up to sign. To avoid notational 
clutter, we use the A-notation to denote both the spanned parallelepipeds and their vol¬ 
ume, computed as the determinant of the matrix with columns given by the factors of 
the wedge product. 

On the one hand, we have 


Uf n = u t - e i 


det u t e 2 



n 

u t A e 2 A • • • A e n = u t A f\ e*. 

i=2 


On the other hand, we first observe that 

det (d^ 0 dT| (t x) ) = detd^detdT| (t x) = detdtf^, 


(7) 


since det d(/?Q = 1 for all t. Finally, we compute det ^d^dTl^ ^ J as the change of volume 
under the action of d</?Qd'F|^ x y 


e t A e 2 A ... A e n 


n n 

d(^o w t A f\ei = u t A f\ei, 

i =2 i =2 


where we have used Eqs. (4) and (5). With Eq. (7), this finishes the proof. 
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B. Proof of det = u • h 


We choose a local unit volume tangent space basis as in Section A. The tangent space 
T(t,x)H i s then spanned by the orthogonal basis (et + fie i, e 2 ,..., e n ). With a = || e* + 
/3ei|| = yjl + fi 2 , we have that ^(et+fie i) Ae 2 A.. .Ae„ is a unit n-volume parallelepiped 

in T (t,x)U- 

Next, we calculate det (dc^d'E) as the change of volume under d^dT) ^ 


det 



( h I (t ,*) ( et + ^ ei )) A AiU d $ d * | (tjX) e* 

d(et + /3ei) Ae 2 A... Ae„ 

^(^(x) - /3ei) A e 2 A ... A e n 
^( u t( x ) — /3ei) • n. 


From the first to the second line, we have used that the denominator is normalized, and 
the different signs are due to Eq. (5). 

On the other hand, let us calculate the n-component of the extended velocity lit = 
et + Ut, i.e., the component normal to the section’s tangent space 

(e t + u t )-h = (e t + w t ) A d ( e * + /3ei) A e 2 A ... A e n 
= d et A (u t - fie i) A e 2 A ... A e n 
= d(w t - fie i) A e 2 A ... A e n 
= d(u t (x) -/?ei) -n. 


From the first to the second line, we have used determinant rules, and from the second 
to the third lines Laplace expansion along the time component. 

Finally, the above calculation is, of course, consistent with the one given in Section A, 
since u • n = u • h whenever n = (^), i.e., when the section does not move. In this case, 
f} = 0 and a = 1. 
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